Quantum Communications Made Easy T 
Deterministic Models of Bosonic Channels 

Graeme Smith and John A. Smolin 
IBM T.J. Watson Research Center, Yorktown Heights, NY 10598, USA 
(Dated: November 9, 2012) 



Information theory establishes the ultimate limits on performance for noisy communication systems 
[T]. An accurate model of a physical communication device must include quantum effects, but 
typically including these makes the theory intractable. As a result communication capacities are 
not known, even for transmission between two users connected by an electromagnetic waveguide 
subject to gaussian noise. Here we present an exactly solvable model of communications with a fully 
quantum electromagnetic field. This allows us to find explicit expressions for all the point-to-point 
capacities of a noisy quantum channel, with implications for quantum key distribution, and fiber 
optical communications. We also develop a theory of quantum communication networks by solving 
some rudimentary quantum networks for broadcasting and multiple access. When possible, we 
compare the predictions of our new model with those of the orthodox quantum gaussian model and 
in all cases we find capacities in agreement to within a constant number of bits. Thus, in the limit of 
high signal to noise ratios our simple model captures the relevant physics of gaussian models while 
remaining amenable to detailed analysis. 

OO 

O , I. INTRODUCTION 



A fundamental property of any communication system is the maximum rate of data transmission possible using the 
best communication schemes. This is called the capacity of a channel. It is usually calculated as a function of noise 
levels and subject to a limited power budget. In 1948, Shannon Q] presented a beautiful theory of information both 
i— — i formulating and solving the capacity problem. For the specific case of a channel with additive white gaussian noise 
his formula can be solved explicitly giving the classical capacity C as a function of the signal to noise ratio (SNR) as 

> 1 

\0 C(SNR) = -Wlog(l + SNR). (1) 

in 2 

This formula guided the development of practical schemes that are now in use, culminating in efficient codes that 
come close to achieving the theoretical limit [5] . 

Noise is not a purely mathematical abstraction, but must arise from some physical process. Such processes are 
properly described, of course, by quantum mechanics, and therefore calculating the true information carrying capacity 
of a channel requires a quantum-mechanical treatment. It is also natural, then, to consider new types of capacities, 
• . such as the capacity of a channel for transmitting quantum states coherently (the quantum capacity), or classical states 
securely (the private capacity). Unfortunately, unlike in classical information theory, for most quantum channels none 
of these capacities are known (see, for for example, O 0] for the quantum capacity and [5] for the classical capacity, 
^ and [5] for the private capacity). 

Whereas much of the existing work on quantum channels has concentrated on abstract finite-dimensional channels, 
here we would like to study the problem in a more realistic setting. Our method is well suited to channels consisting of 
gaussian noise in bosonic electromagnetic modes, though it is substantially more general. We will be able to calculate 
classical, quantum, and private capacities for a wide range of realistic channels. 

Because quantum information cannot be cloned [7; , knowledge gained by the environment about a signal is necessar- 
ily detrimental to quantum transmission. This need to consider what information is transmitted to the environment 
as well as what goes to the intended receiver puts an analysis of quantum capacity on par with the study of the clas- 
sical multi-user broadcast channel (which is also notoriously difficult to analyze [5]). To make the problem tractable 
we, following the work of [5J, substitute a discretized and deterministic model for the actual channel. In the limit 
of high SNR this model will capture the important features of the real channel, and allow us to calculate capacities 
to within a small number of bits. The general approach can be thought of as discretizing the continuous system 
under consideration in a very simple way and then truncating signals smaller than the noise power. The result is a 
deterministic channel model that is easy to analyze. 
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II. ADDITIVE GAUSSIAN CHANNELS 

States of electromagnetic modes are described by their quadratures P, Q and are called gaussian when they can 
be completely characterized by a matrix 7 containing the covariances of these quadratures |10j . Such states can be 
visualized by their Wigner functions — quasi-probability distributions depicting the state's location in phase space. 
Gaussian states have Wigner functions which are ellipsoids with gaussian profiles. Fig. [l|a) shows the Wigner function 
for a gaussian state with covariances ap, o~q. The minimum uncertainty state has <jp<jq = 1/2 (in units of K). Such 
states are always pure. Mixed states have (Jpoq > 1/2 and can be thought of as mixtures of pure states. 

We will replace the quadratures P,Q with a simpler discretized model described in Fig. [ljb). We call this the 
discrete quadrature (DQ) model. Roughly speaking, states in phase space are more distinguishable when their Wigner 
functions are less overlapping. To reflect this property, we divide phase space up into a lattice of nonoverlapping 
rectangles that we take to be perfectly distinguishable. The smallest physical state has area 1/2. The Wigner function 
is thus replaced with several rectangles of area 1/2 which tile the region of phase space where it is nonnegligable. The 
associated state in our model is a uniform mixture of these distinguishable rectangles. 

The model of the action of a channel model will also be discrete: Since all the output states are perfectly distinguish- 
able rectangles, every pair of input states will either be mapped to distinguishable outputs, or to the same output. 
Calculating the achieved communication rate is then a simple matter of counting the distinguishable output states. A 
full calculation of the capacity will also include a maximization over modulation schemes, i.e. the set of input states 
employed. Note that not all modulation schemes are physical. The choice must obey the following constraints: 

• K a Q + °p) — W ( nmte power) 

• APAQ > 1/2 (quantum uncertainty) 

• AP < ap, AQ < aQ (common sense) 
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FIG. 1: Phase-space representation of states (Wigner Functions): (a) Wigner function of a gaussian mixed state with variances 
aQ and ap. (b) The discrete quadrature model. The mixed state has been approximately decomposed into nonoverlapping 
rectangles with width AQ and height AP. These rectangles are the discrete states of the model which replaces original physical 
system. We imagine them as approximations to the pure squeezed state shown at point A. 

The evolution of bosonic states under the action of channels is described by the evolution of the system's quadratures. 
For gaussian noise, which arises from quadratic Hamiltonian interaction with gaussian environment modes, the allowed 
evolutions take a particularly simple form, being completely described by a symplectic matrix |llj . For example, 
P — » vAP + y/1 — X r, Q — » vXQ + \/l — A s, results from interaction with an environment mode described by r, s 
(see Fig. [2]). The additive gaussian noise channel arises when the environment begins in a gaussian state, and the 
thermal noise channel, when it is initially in a thermal state. 
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FIG. 2: Additive gaussian noise noise channel. An input state in beam A is combined with an environmental input pE on a 
beamsplitter with transmitivity A. One output beam is discarded (remains available only to the environment) and the other is 
the channel output B. When the environment's input pE is a vacuum state, this gives the attenuation channel. When pE is a 
thermal state, it is the thermal noise channel. 

A. Classical Capacity 

1. General Additive Gaussian Noise 

We now apply our discretization procedure to model classical communication over an additive gaussian noise 
channel. Fig. [3] illustrates the following analysis. We suppose that the input signal has some average power constraint 
W > \{&p + Oq). Furthermore, we fix a modulation scheme, to be optimized over later, which amounts to deciding 
the shape of the rectangles in the discretization shown in Fig. PIT). Given the variables <7p,(jQ and rectangle shape 
AP, AQ with APAQ = 1/2, our input space has (Jp<7q/(APAQ) distinguishable states (Fig. [3^i). We must now 
determine how many distinguishable outputs these get mapped to. Attenuation by A maps the entire input space to 
a VXap x %/A(7q rectangle of y/XAP x \/AAQ tiles (Fig. K3p). When a noise of typical size \fl — Xa r is added to the 
P quadrature, tiles closer than this are taken to be confusable, and similarly for Q, which leads us to "meta"-tiles of 
dimension 

m&x(VXAP, Vl - Aoy) x max(v / AAg, VI - Xa s ) (Fig. ^p). (2) 

If it happens that this rectangle's area is less than 1/2 the tiles in Fig. [3fc are smaller than the minimum tile size 
allowed by uncertainty. We must then chose a tile shape, AN x AM with AN AM — 1/2 satisfying uncertainty. We 
thus find a final tile dimension of AP x AQ with 

AP = max( ^A AP, Vl - Xa r , AN) (3) 
AQ = max(^AAQ, \fl - Xa s ,AM) (Fig.|§l). (4) 

This gives a total number of distinguishable output states of Ao\p<7q/(APA(5) and a classical capacity C of 

C = maxlogf^^ (5) 
B \APAQj y ' 

where the maximization is over all of the constrained variables: \{cfp + Cq) < W, APAQ > ^, AP < <jp, AQ < 
ctq, AN AM > We now evaluate the preceding formula for some important special cases. 

2. Example: Attenuation 

The attenuation channel is an additive channel with transmissivity A and pure environment in a vacuum state with 
a r = at = -4? . To good approximation, this channel describes the propagation of the signal through lossy optical 



4 



a 



a P 



AQ 



AP 



Attenuation 



XAQ 

Taap 



(b) 



(d) 



AQ = max (vAAQ, s/l^Xa s , AN\ 

AP = max ( v/A AP, VT^Vr , AM 



Noise 



max ( x/AAQ, \/l — Aer s 



Heisenberg 



\/A, 



dp 



VAAP yT^Aov) 



\/A<xq 



FIG. 3: The DQ model of the action of a channel, (a) The input space is divided up into tiles of AQ x AP. (b) Attenuation 
has shrunk the tiles to \T\AQ x \/AAP. (c) Noise has redenered some tiles confusable, so the effective tile size is increased, 
(d) If any tiles are smaller than the smallest tile size allowed by uncertainty, the tiles are enlarged again. 



fiber. Letting AP = AQ 



-j= and AM = AN = minimizes APAQ = | while satisfying the constraints. This 

leads to a capacity of log|~2V^A~|. For this channel, we actually know the classical capacity exactly [12]. Indeed, the 
true capacity g (A(W — o)) / l n 2, where g(x) = (x + 1) ln(a; + 1) — x hix, differs from our estimate by no more than 
1.4 bits (see Fig. 



3. Example: Classical Noise 

When a channel applies gaussian-distributed kicks in phase space, it is called a classical noise channel. This arises as 
a limiting case of the thermal noise channel shown in Fig. [2]with A — > 1 and a r = a s — > oo while keeping a r s/l — A = /i. 
We would like to evaluate Eq. ^ for this channel, where we have AP = max(AP, /x) and AQ = max(AQ, jj,). 

First we consider the case p? > 1/2. First note that APAQ > p? . Choosing AP = AQ = l/v2 achieves this 
lower bound. Note also that max H/> 1^2 &pQq — W is achieved for up = <jq = VW. Thus by choosing 

AP = AQ = 1/V2 and up = oq = \/W we simultaneously maximize opOQ and minimize APAQ while satisfying 
all the constraints in Eq. We thus find a capacity of log(W/ fj 2 ). In a similar way, when /1 2 < 1/2, APAQ is 
minimized by AP = AQ = l/v2 so that by the same argument the capacity is now log(2M^). A lower bound for the 
capacity of the classical noise channel with noise power \j? is (g(W - 1/2 + p 2 ) - g{p 2 )) / In 2. This is achieved by 
classical displacements of vaccuum states. This bound is within 1.45 bits of the true capacity [T3] and agrees with 
the DQ model to within 1 bit for all . 
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FIG. 4: The classical capacity of the attenuation channel, (a) For A = 1/2, the actual capacity g((W — l/2)/2)/ln2 and the 
capacity calculated using the discrete quadrature model log W. (b) The difference between the actual and discrete quadrature 
capacity as a function of both A and power. Note that the difference is never greater than 1.4 bits. 




(a) (b) 



FIG. 5: Dephasing channel, (a) Discrete quadrature model. The max[AQ, fj] in the denomiator of the Q quadrature output 
takes into account that you cannot distinguish states smaller than either the size of the input or the size of the applied classical 
noise. The 'ceiling' brackets inside the logarithms avoid them blowing up if the noise is so high that there are fewer than one 
distinguishable states by this calculation, (b) Classical capacity of the dephasing channel with fi — 1. The capacity is achieved 
with AP = l/(2\/W) and AQ — VW and op — oq = VW. The squeezing and therefore the capacity is limited only by the 
input power. The lower line in the plot is the prediction from the discrete quadrature model, the upper is an exact calcuation 
of the Holevo quantity given a signal ensemble of maximally-squeezed pure states subject to the power constraint, which is 

an achievable rate \ = 2W/ ^ 1+M 2 /( 2W > 1 ^ _ g(^A +g M 2W ) / m 2. Our model suggests this is not far from the actual 
capacity. 



4- Example: Dephasing 

Another interesting special case is the classical dephasing channel, which adds classical noise of power /i 2 to the 
Q quadrature while leaving P untouched. Starting from the quantum channel P — > V XP + s/i — Xr and Q — > 
VXQ + \/l - As, we take the limit A — >• 1 with (1 — X)a 2 s = fi 2 , (1 — A)er 2 -> 0. So we have the channel P -> P and 
Q — > Q + £ where £ is a classical guassian random variable with variance /i 2 . 

The discrete quadrature model of this channel is considered in Fig. [5] , where we have obtained a capacity estimate 
of log2VK, independent of noise level. The modulation scheme to achieve this involves highly squeezed signal states. 
By choosing states very narrow in P and broad in Q, we can pack nearly all of the signal into the noiseless quadrature. 
Thus, for classical transmission, such a channel achieves rates just as high as in the noiseless channel, albeit with more 
difficult modulation. This prediction is borne out by computing the Holevo information on the suggested squeezed 
ensemble, which an gives achievable rate that differs from log 2W by no more than a bit [TTJ Q3] . 
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FIG. 6: (a) Decomposition of thermal noise channel with average photon number Ne into a pure attenuation channel with 
A = \/G and a gain channel with gain G. (b) DQ model of this decomposition. 

B. Quantum Capacity 

As mentioned above, evaluating the quantum capacity Q requires assessing not only signals sent from sender 
to receiver, but also how information is leaked from sender to the channel's environment. While this makes our 
calculations more complex, we will nevertheless find simple and reliable estimates for quantum capacities. The 
quantum capacity will be computed as the number of distinguishable states communcated to the channel's output 
about which the environment knows nothing at all. This complete lack of knowledge enables one to communicate 
quantum superpositions of these distinguishable states, therefore our model imagines them to define the basis states 
of a Hilbert space that will be successfully transmitted. Note that this is exactly the definition of the private capacity 
|15j . and therefore our model will not be able to separate quantum and private capacities. 

1. Example: Thermal Noise 

We now turn our attention to the thermal noise channel, which maps 
where r, s are the quadratures of a thermal state with average photon number Ne- Such a channel can be decomposed 
as a composition of a pure loss channel with an ideal amplifier as shown in Fig. [6^. This decomposition allows us to 
express the environment's state as 

yr^i^Q)-^!,^) (6) 

VG(r 2 ,s 2 ) + VG^l(P,-Q) (7) 

where G = (1 — X)Ne + 1 and A = A/G and r\, si and r 2 , s 2 are the quadratures of two independent vacuum states. 

Following our prescription, an input with power W and variances (jp + Oq < W we have X<jp<jq/ APAQ distin- 
guishable states at the output, or 

Output = \og(Xa P a Q /APAQ) (8) 

bits. The highest order bits of the input get mapped to log(Aerp<7Q/APA(3) states at the output. Similarly, there 
are two environmental modes. Similarly, according to Eq. |6]), the first environmental mode gets 

Si = log ((1 - X)a P a Q /AP 1 AQi) - log (2(1 - X/G)a P a Q ) (9) 

bits where AP\ = max( y/l — XAP, AM\ , \fXa ri ) and AQi = max(\/ 1 — AAQ, AiVj., \/~Xa Sl ). We always have 
AP 1 AQ 1 > 1/2, and AP 1 AQ 1 = 1/2 can be achieved for appropriate AM b AN X . 

For the second environmental mode, we have S 2 — log ((G — l)ap(jQ/AP 2 AQ 2 j = log(2(G— 1)ctp(Tq/G) where 
AP 2 = max(VG - 1AP, AM 2 , VGa r J and AQ 2 = max(VG - 1AQ, AiV 2 , VG<t S2 ) and AP 2 AQ 2 = G/2 is achievable. 
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FIG. 7: Simultaneous classical and quantum communication over a gaussian attenuation channel, (a) DQ Channel model. 
The higherst-order (blue) bits are seen by the environment and therefore can never be used for quantum communication. They 
can always be used for classical communication, even when the lower order bits are being used for quantum communication, 
(b) Rate region compared to the time-sharing tradeoff. 



Note that the high-order bits transmitted to both mode 1 and mode 2 are identical, but mode 1 gets more of them 
so that the total number of bits leaked to the environment is S\ . The quantum capacity is therefore 

S ou tput ~Si = log(Xa P a Q /APAQ) - log (2(1 - \/G)a P a Q ) (10) 
= logA-log(l-A/G)-log2APAQ (11) 

- log A - log(l - (1 _ A) A ^ + 1 ) - log 2APAQ (12) 



We now evaluate ( 12 ) for the attenuation channel (Ne =0): As we showed in our discussion of the classical capacity 
of the attenuation channel, by appropriate choices of AP, AQ and AM, A7V, we can achieve APAQ = 1/2. Then the 
capacity is log A — log(l — A) for A > 1/2 and unlimited input power. This is exactly the quantum capacity of the 
gaussian attenuation channel |16j . In fact our model offers slightly more information. Throughout in order for our 



model to make sense, the estimates of the number of levels transmitted must be integers. To ensure this in Eq. ( 10 ) 
we need to have power W > 1/A, which suggests the minimum power necessary to achieve capacity. 



III. SIMULTANEOUS QUANTUM AND CLASSICAL COMMUNICATION 

Our model also allows for a simple understanding of the recently discovered [17] tradeoff between classical and 
quantum communication over the same optical channel. Switching between the two different optimal communication 
strategies for the two types of information (time sharing), gives a linear tradeoff between the two communication 
rates. In |17j . it was shown that rate region for an attenuation channel could be substantially larger than this naive 
strategy would suggest. A glance at Figure [7] makes it clear that while sending quantum information at capacity, it 
is possible to simultaneously send some classical information with no degradation of the quantum transmission. Only 
after sending classical information at a rate greater than C — Q docs a linear tradeoff emerge. 
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(a) (b) 
FIG. 8: Multi-user Channels: (a) Multiple-access channel, (b) Broadcast channel. 



IV. MULTI-USER COMMUNICATION 

A useful and natural generalization of the channel capacity problem is when there are multiple senders and/or 
receivers. Our simple model is well adapted to extracting useful answers in this setting which is typically highly 
intractable, even in the classical setting [151 US]- Below we consider two multi-user channels: The multiple access 
channel with two senders and one receiver and the broadcast channel with one sender and two receivers (see Fig. |8| . 

A simple multiple access channel is shown in Fig. [8^,. Two senders A and B try to transmit information to a 
reciever C . The performance of such a channel is described by a rate region rather than a single capacity. The 
discrete quadrature model of this channel is shown in Fig. [9£t and the predicted rate region is compared to existing 
bounds in Fig. [9]3. Similarly, the DQ model and rate region for a simple broadcast channel, where one sender, A, 
tries to send information to two receivers B and C, is shown in Fig. [lO} 




Ra 

(a) (b) 

FIG. 9: Simple multiple access channel, (a) Discrete quadrature model for the multiple access channel from Fig[8ji with XWa > 
(1 — \)Wb- Discretizing the input quadratures as usual gives \ogap A <jQ A / APaAQa for sender A and log <Jp b oq b / APp AQb 
for sender B. The channel maps input quadratures Pa,Qa and Pb,Qb to outputs V\Pa,Qa + \/l — \Pb ,Qb- The highest 
order bits from A remain distinguishable at the output regardless of what B sends. The lower order bits are confusable; one of 
the senders can get information through to C reliably so long as the other sender holds the corresponding bits fixed, (b) The 
classical-capacity rate region for the multiple access channel. The outer blue region is the upper bound from [20]. The black 
lines indicate their achievable rate region with coherent states and the achievable region for our model with APa = AQa, 
APb = AQb- Our region is lower by just a fraction of a bit. 



V. DISCUSSION 

Our approach allows us to solve a number of vexing questions that are intractable in the orthodox model. For 
example, we can exactly solve models for single-user communication with thermal noise, as well as multi-user networks 
including broadcast and multiple- access channels. Perhaps, though, our model is too crude to capture the relevant 
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FIG. 10: Broadcast channel, (a) The discrete quadrature model for the broadcast channel of |8p. We have taken AP = AQ = 
1/V2 which achieves optimality. The rate region is then given by R c < log 2(1 — X)W and Rb + Rc < log 2AM 7 . (b) The rate 
region is plotted for W — 50.5, A = .8 and compared to the lower bound from 21 which they conjecture to be the actual 
capacity. 
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behavior of quantum communication systems. We argue, to the contrary, that we capture the relevant physics by 
comparing the few solvable gaussian examples to our predictions where we find good agreement to within one or two 
bits. Furthermore, for examples where only lower bounds are available, we find good agreement with these, thus 
predicting that typically known lower bounds are equal, or at least close, to the ultimate capacities. 

Our model also explains some previously known but counter-intuitive facts, rendering them almost obvious. For 
instance, we can explain why while classical capacity rises without bound as power increases, the quantum capacity 
saturates: Increasing power enhances transmission to both receiver and eavesdropper in equal measure. Put simply, 
if you're trying to transmit privately, shouting your secrets doesn't help. We have explained the results of |17j . that 
time-sharing is not optimal for the classical-quantum tradeoff when trying to simultaneously transmit some of each 
type of information. 

Finally, we can make some predictions. Within our deterministic model, entanglement between channel uses and 
other quantum tricks don't appear to be useful. In particular, (1) Privacy and coherence are equivalent in our 
model. Therefore we predict the private and quantum capacities will always be nearly equal in gaussian channels, 
even though they can be very different in general [22 j . (2) Two-way communication doesn't help much for gaussian 
channels, again counter to the general case [23] . (3) Gaussian channels have a single-letter capacity forumla to within 
a small number of bits (compare to [231 [35]). (4) Capacities of gaussian channels are nearly additive (unlike the 
extreme superadditivity in ;24, 125]). 

The "nearly" can, however, obscure some interesting effects. We know, for example, that gaussian channels can 
display superactivation, that is there exist pairs of gaussian channels each with zero capacity that can neverthelesss 
be used together to achieve positive capacity The resolution is that the resulting capacities are very small (the 
rate achieved with the joint channel in [26] is only 0.06 bits). Our predictions are only meant to be with a few bits 
of the "correct" value so there is no contradiction, and our results should be asymptotically correct at high SNR. 
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